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Some numerical characteristics of bipartite graphs in relation to the 
problem of finding all disjoint pairs of S-permutation matrices in the gen- 
eral X case are discussed in this paper. All bipartite graphs of the 
Q ! type g = {RgU Cg,Eg), where \Rg\ = \Cg\ = 2 or \Rg\ = \Cg\ = 3 are 

provided. The cardinality of the sets of mutually disjoint S-permutation 
matrices in both the 4x4 and 9x9 cases are calculated. 
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^ ; 1 Introduction 

' Let m be a positive integer. By [m] we denote the set 

p. H = {l,2,...,m}, 

. We let Sra denote the symmetric group of order ra i.e., the group of all one- 

' to-one mappings of the set [to] to itself. If x G [to,], p G 5m, then the image of 

the element x in the mapping p we will denote by p{x). 
A bipartite graph is an ordered triple 



9 = {Rg,Cg,Eg), 

where Rg and Cg are non-empty sets such that Rg CiCg = 0. The elements of 
Rg U Cg will be called vertices. The set of edges is Eg C Rg x Cg = { (r, c) | r e 
-Rg, c e Cg}. Multiple edges are not allowed in our considerations. 

The subject of the present work is bipartite graphs considered up to isomor- 
phism. 

We refer to [3] or [6] for more details on graph theory. 

Let n and k be two nonnegative integers and let < fc < n^. We denote by 
<Sn,k the set of all bipartite graphs of the type g = {Rg,Cg, Eg), considered up 
to isomorphism, such that |i?g| = |Cg| = n and |i?g| = k. 

Let Pij, i i, j < ri, be square nxn matrices, whose entries are elements 
{1,2,..., n^}. The x matrix 



1 



p = 



Pll Pl2 
P2I P22 



Pin 
P2n 

P 



is called a Sudoku matrix, if every row, every column and every submatrix 
Pij, 1 < i, j < n comprise a permutation of the elements of set [n^], i.e., every 



number s e {1, 2, . . . 
Pij. Submatrices Pij 



n?} is found just once in each row, column, and submatrix 
are called blocks of P. 
Sudoku is a very popular game and Sudoku matrices are special cases of 
Latin squares in the class of gerechte designs [I] . 

A matrix is called binary if all of its elements are equal to or 1. A square 
binary matrix is called permutation matrix, if in every row and every column 
there is just one 1. 

Let us denote by E„2 the set of all x permutation matrices of the 
following type: 



A = 



An 
A21 



A12 
A22 



Anl A„ 



Aln 
A2n 

A 



where for every s,t g {1,2,. Ast is a square n x n binary submatrix 
(block) with only one element equal to 1. 

The elements of E„2 will be called S-permutation matrices. 

Two E„2 matrices A — {oij) and B = (bij), 1 < i,j < will be called 
disjoint, if there arc not elements Oij and bij with the same indices such that 



Oiij — bij — 1. 



The concept of S-permutation matrix was introduced by Geir Dahl 2J in 
relation to the popular Sudoku puzzle. 

Obviously, a square xn^ matrix P with entries from [v?] = {1, 2, . . . , n^} is 
a Sudoku matrix if and only if there are E„2 matrices Ai, A2, 
disjoint, such that P can be written in the following way: 



, Ari 



pairwise 



P=l- Ai + 2- A2 



n" ■ A„2 



(1) 



In [S] Roberto Fontana offers an algorithm which returns a random family of 
X mutually disjoint S-permutation matrices, where n = 2, 3. For rt = 3, he 
ran the algorithm 1000 times and found 105 different families of nine mutually 
disjoint S-permutation matrices. Then, applying (P), he decided that there are 
at least 9! • 105 = 38 102 400 Sudoku matrices. This number is very small 
compared with the exact number of 9 x 9 Sudoku matrices. In [3] it was shown 
that there are exactly 

9! • 72^ •2^-27 704 267 971 = 6 670 903 752 021 072 936 960 

number of 9 x 9 Sudoku matrices. 

To evaluate the effectiveness of Fontana's algorithm, it is necessary to cal- 
culate the probability of two randomly generated matrices being disjoint. As is 
proved in [2^ , the number of S-permutation matrices is equal to 



(n!) 



2n 



2 



Thus the question of finding a formula for counting disjoint pairs of S- 
permutation matrices naturally arises. Such a formula is introduced and verified 
in [5]. In this paper, we demonstrate this formula to compute the number of 
disjoint pairs of S-permutation matrices in both the 4x4 and 9x9 cases. 

2 A formula for counting disjoint pairs of S- 
permutation matrices 

Let g = {Rg, Cg, Eg) G &n,k for some natural numbers n and k and let v E Vg ^ 

Rg U Cg. 

By N{v) we denote the set of all vertices oiVg, adjacent with v, i.e., u e N{v) 
if and only if there is an edge in Eg connecting u and v. If w is an isolated 
vertex (i.e., there is no edge, incident with v), then by definition N(v) = and 
degree(u) = \N{v)\ =0. If v £ Rg, then obviously N{v) C Cg, and if u G Cg, 
then N{v) C Rg. 

Let g = {Rg,Cg, Eg) e 0n,k and let u,v e Vg — Rg U Cg. We will say that 
u and V are equivalent and we will write u ^ v if N{u) — N(y). If u and v are 
isolated, then by definition u ^ v ii and only if u and v belong simultaneously to 
Rg, or Cg. The above introduced relation is obviously an equivalence relation. 

By Vg^^ we denote the obtained factor-set (the set of the equivalence classes) 
according to relation ^ and let 

K,/^ = {Ai,A2,...,AJ, 

where A, C Rg, or A, C Cg, i = 1,2, ... s, 2 < s < 2n. We put 

Si = \Ai\, l<6i<n, j = l,2,...,s 

and for every g £ &n,k we define multi-set (set with repetition) 

[g] = {Si,S2, ...Ss}, 

where 61,62, ... ,5s are natural numbers, obtained by the above described way. 

li Zi Z2 ... Zn is a permutation of the elements of the set [n] = {1, 2, . . . , n} 
and we shortly denote p this permutation, then in this case we denote by p{i) 
the i-th element of this permutation, i.e., p{i) ^ Zi, i = 1,2, ... , n. 

The following theorem is proved in [5]: 

Theorem 1 ^81 Let n > 2 be a positive integer. Then the number Z?„2 0/ all 
disjoint ordered pairs of matrices in I]„2 is equal to 

= (n!)^" + (n!)^("+i) ^(-1)'= ^ ^^^.-c 

fc=i se©„,fc 11 '5' 

"5e[g] 

The number d„2 of all non-ordered pairs of disjoint matrices in S„2 is equal 
to ^ 

□ 
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The proof of Theorem [T] is described in detail in [5] and here we wih miss it. 

In order to apply Theorem [T] it is necessary to describe all bipartite graphs 
up to isomorphism g ~ {Rg,Cg, Eg), where \Rg \ — \Cg \ — n. 

Let n and k are positive integers and let g G (Sn,fc- We examine the ordered 
{n + l)-tuple 

^{g)^{M9),M9),---,M9)), (4) 
where ipiig), i = 0, 1, . . . , n is equal to the number of vertices of g incident with 

n 

exactly i number of edges. It is obvious that iipiig) — 2fc is true for all 
9 G &n,k- Then formula ^ can be presented 

A.^ = (r.!)4" + (n!f(«+i)^(_lf- ^ . 

k=l gG®„.fc 11 

sela] 

Since (n - n)! = 0! = 1 and [n - (n - 1)]! = l! = 1, then 

n-2 



i'i(g) 



D„. = (n!)4" + (n!f("+i)^(-l)'= ^ . (5) 

fc=i gee„.fc 11 
ma] 

Consequently, to apply formula © for each bipartite graph g £ &n.k and for 
the set &n.k of bipartite graphs, it is necessary to obtain the following numerical 
characteristics: 

n-2 

ni(..-i)!i*''" 

selg] 

and 

0(n,fc)= ^ uig) (7) 

Using the numerical characteristics ^ and ([7]), we obtain the following 
variety of Theorem [TJ 

Theorem 2 

A.^ - (n!)4" + (n!)2("+i) fc), (8) 

fe=i 

where 9{n, k) is described using formulas ^ and @). 

□ 
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3 Demonstrations in applying of Theorem [2] 



3.1 Counting the number D4 of all ordered pairs of disjoint 
S-permutation matrices for n = 2 

3.1.1 k = l 

h\ n — 2 and k — 1, ©2.1 consists of a single graph gi shown in Figure [T] 




Figure 1: n = 2, k = 1 
For graph gi G ©2.1 we have: 

N = {1,1, 1,1} 

= (V'o(5i),^i(5i),^2(<?i)) = (2,2,0) 

Then we get: 



and therefore 



, , [(2-0)!]^ ^ 

0(2,1)= ^ ^(5) = 4. 
see2,i 



(9) 



3.1.2 k^2 

The set ©2,2 consists of three graphs 172, 93 and f/4 depicted in Figure [2] 






Figure 2: n = 2, fc = 2 
For graph 172 G ©2,2 we have: 

[52] = {1,1, 1,1} 
^(52) = {M92),M92),M92)) = (0,4,0) 

[(2-0)!]" 
^^^^^ " 1!1!1!1! ^ ^ 
For graphs 53 G ©2,2 and 174 G ©2,2 we have: 

[53] = M = {2,1,1} 
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*(33) = *(54) = (1,2,1) 
w(33) = c.(g4) = ifcS^ = 1 



2! 1! 1! 



Then for the set ©2,2 we get: 



e{2,2)^ J2 0.(3) = 1 + 1 + 1 = 3. 

S6«2,2 

3.1.3 fc = 3 

In n = 2 and fc = 3, ©2.3 consists of a single graph shown in Figure |31 



(10) 



95 



Figure 3: n = 2, k — 3 
For graph 35 G ©2,3 we have: 

[55] = {1,1, 1,1} 

*(55) = (V'0(55),^l(55),V'2(g5)) = (0,2,2) 

Then we get: 



and therefore 



3.1.4 fc = 4 



[(2-0)!]° 
^^^'^ = 1!1!1!1! = ^ 

0(2,3) = ^ c^(5) = l. 

g6e2.3 



(11) 



When n = 2 and fc = 4 there is only one graph and this is the complete bipartite 
graph which is shown in Figured 




Figure 4: 71 = 2, fc = 4 
For graph G ©2,4 we have: 

[56] = {2, 2} 

*(56) = {Mg6),M96),M96)) = (0,0,4) 
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Then we get: 



[(2-0)!f _ 1 



2! 2! 4 

and therefore 

0(2,4)= io{g) = ^. (12) 

se«2,i 

Having in mind the formulas dH]), ®, (UHl), (HU and ([H]) for the number D4 
of ah ordered pairs disjoint S-permutation matrices in n = 2 we finally get: 

Di = (2!)8 + (2!)6 [-9(2, 1) + 0(2, 2) - 0(2, 3) + 0(2, 4)] = (13) 

1^ 



256 + 64 -4 + 3 -1 



144. 



The number of all non-ordered pairs disjoint matrices from E4 is equal to 



di ^ -Da = 72. 
2 



(14) 



3.2 Counting the number Z^g of all ordered pairs of disjoint 
S-permutation matrices for n = 3 

3.2.1 k = l 

Graph 57, which is displayed in Figure [5] is the only bipartite graph belonging 
to the set ©3 1. 



97 
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1^ J 


Rg7 





Figure 5: n — k — 1 
For graph 57 G ©3.1 we have: 

M = {1,1,2,2} 

^(57) = (V'o(<?7),V'l(<?7),^2(57),V'3(57),V'4(58)) = (4,2,0,0) 

Then we get: 



^^(57) 



[(3-0)!]^[(3- l)f _ 64-22 



1! 1! 2! 2! 



1 • 1 • 2 • 2 



= 1296 



and therefore 



0(3,1)= ^(3) = 1296. 

se©3,i 



(15) 
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3.2.2 k = 2 



In this case 03^2 — {58:ff9;3io}- The graphs g^, gg and giQ are shown in Figure 

El 
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^910 



Figure 6: n = 3, fc = 2 
For graph gs G ©3,2 we have: 

M = {1,1,1,1,1,1} 

*(58) = {M9s),M9s),M98),M98),M98)) = (2,4,0,0) 

l! l! l! l! 1! 1! 

For graphs 59 G ©3,2 and 510 £ ©3,2 we have: 

[99] - [510] -{1,1,2,2} 
vl/(59) = *(ffio)-(3,2,l,0) 



^^(59) = w(5io) 



[(3-0)!]^[(3-l)!]' _ 63-22 



1! 1! 2! 2! 



1-1-2-2 



= 216 



Then for the set ©3^2 we get: 



61(3,2)= w(.g) = 576 + 216 + 216 = 1008. 

ge@3.2 



(16) 



3.2.3 k = 3 

When n = 3 and fc = 3 the set ©3,3 = {511, 512, ffi3, 5i4, 315, ffie} consists of six 
bipartite graphs, which are shown in Figure [T] 
For graph gu e ©3,3 we have: 

[gn] ={1,1,1,1,1,1} 
vl'(5n) = (0,6,0,0) 

[(3-0)!]°[(3-l)f _ 
'^(511)- -6 -2-64 
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Figure 7: n = 3, fc = 3 

For graphs 512, 5i3 € ©3,3 we have: 

M = [513] = {1,1, 1,1, 2} 

f(5l2) = *(5l3) = (1,4,1,0) 

/ ^ [(3-0)!]M(3- l)f 61-24 
For graph gii4 e ©3,3 we have: 

[514] = {1,1, 1,1, 1,1} 

*(5i4) = (2,2,2,0) 



For graphs 515,516 G ©3,3 we have: 

[515] = [516] = {1,2,3} 

*(5i5) = *(5i6) = (2,3,0,l) 



1!2!3! 



2-6 



Then for the set ©3.3 wc get: 



6»(3,3)= ^ w(5) = 64 + 48 + 48 + 144 + 24 + 24 = 352. (17) 
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3.2.4 fc = 4 



When n = 3 and /c = 4 the set 63^ = {517,518,919,520,521,522,523} consists of 
seven bipartite graphs, which are shown in Figure [5] 




718 
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Figure 8: n = 3, fc = 4 
For graph 317 e ©3^4 we have: 

[517] = {1,1, 2, 2} 

4'(5i7) = (2,0,4,0) 



[(3-0)f [(3-1)!]" 6^.20 

W(5l7) = = ^r, = 9 



1! 1! 2! 2! 

For graph gig, e ©3^4 we have: 

M = {1,1,2,2} 



22 



uj{g 



18) 



vl/(5i8) = (0,4,2,0) 
[(3-0)!]" [(3-1)!]^ 6" -2^ 



1! 1! 2! 2! 



22 



For graph gig e ©3^4 we have: 



10 



[ffig] ={1,1,1,1,1,1} 



vl/(5i9) = (0,4,2,0) 

^(...)-'^^T??7y^'-60.2^ = 16 
1! l! l! l! l! l! 

For graphs 520 G ©3,4 and 521 £ ©3,4 we have: 

[.920] = [.921] = {1,1,1,1,1,1} 

^(520) = *(52i) = (1,2,3,0) 



[( 3-0)!]M(3-l )!] 
1! 1! 1! 1! 1! 1! 
For graphs 522 £ ©3,4 and 523 £ ©3,4 we have: 



[522] = [<?23] = {1,1,1,1,2} 
*(.922) = *(523) = (1,3,1,1) 
, s [(3-0)!]M(3-l)!]^ 61-23 

W 922 = 923 = -^^^ / \ ' ^ = = 24 



Then we get: 

6»(3,4)= ^ a;(5) = 9 + 4+ 16 + 24 + 24 + 24 + 24 = 125. (18) 

96®3,4 

3.2.5 fc = 5 

When n = 3 and k — 5 the set ©3^5 consists of seven bipartite graphs (724 + ff30, 
which are shown in Figure [SI 
For graph (724 £ ©3,5 we have: 



[524] = {1,1,2,2} 
vI/(524) = (0,4,0,2) 



[(3-0)!]°[(3-l)!f _6°.24_ 
'^(ffis) 111! 2! 2! - ^ 

For graph 325 £ ©3,5 we have: 

[325] = {1,1,2,2} 

vI/(525) = (0,2,4,0) 
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Figure 9: n = 3, A; = 5 



^(518) 



[(3-0)!]° [(3-1)!]' 60-22 



1! 1! 2! 2! 



22 



For graph g2% G ©3,5 we have: 

[526] = {1,1, 1,1, 1,1} 

*(ff26) = (0,2,4,0) 



W(5f26) = 



[(3-0)!]"[(3-l)!]^ _ 2 
1! 1! 1! 1! 1! 1! 



For graphs 527 G 253,5 and 528 € ©3,5 we have: 

[527] = [528] = {1,1, 1,1, 1,1} 
*(527) = *(528) = (0,3,2,1) 

u,(n )-u,(n ^ [(3-0)!]°[(3-l)!]^ _.o 

0J{927) - W(528) - 1, 1, 111! 111! - 6 

For graphs 529 S 63,5 and 530 e 63,5 we have: 
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[529] = [<?30] = {1,1,1,1,2} 



^-(329) = *(530) = (1,1,3,1) 



[(3-0)!]M(3-l)!]^ 61-21 
w(g29) = w(g3o) = TT-rr-rr-Ti-r. = — - — = 6 



Then we get: 



3.2.6 fc = 6 



1! 1! 1! 1! 2! 



0(3,5)= ^(5)-4 + l + 4 



6 + 6 = 37. 



(19) 



When n — 3 and fc = 6 the set 63,6 = {ffsi, 532, ff33, 534, ffss, ffae} consists of six 
bipartite graphs, which are shown in Figure [TUl 








Figure 10: n — 3, k — 6 
For graph 531 e ©3^6 we have: 

[331] = {1,1,1,1,1,1} 




^^(.931 



*(g3i) = (0,0,6,0) 



[(3-0)!]° [(3-1)!]" 



1! 1! 1! 1! 1! 1! 

For graphs 532 £ ©3,6 and 533 £ ©3,6 we have: 

[532] = [333] = {1,1,1,1,2} 



*(532) = *(533) = (0,1,4,1) 
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[(3-0)!]° [(3-1)!]^ 6«.2i 



For graph 334 e ©3^6 we have: 

[534] ={1,1,1,1,1,1} 



vl/(534) = (0,2,2,2) 



34; 



[(3-0)!]^ [(3-1)!]' 6" -22 



1!1!1!1!1!1! 1 

For graphs 535 G 63,6 and 535 6 ©3,6 we have: 

[535]- [936] -{1,2,3} 
^'(535) - *(536) - (1,0,3,2) 

[(3-0)!]M(3-l)!]° _ 61 -20 _ 1 



^(335) = W(g36) = 

Then for the set ©3.6 we get: 



1! 2! 3! 



61(3,6)= J2 ^^(.9) = 1 + 1 + 1 + 4 

Se03,6 



2-6 

1 1 

2 + r 



(20) 



3.2.7 fc = 7 



When n — 3 and k = 7 the set ©3,7 = {.937,338,339} consists of three bipartite 
graphs, which are shown in Figure [TT] 






Figure 11: n — 3, k ~ 7 
For graph 537 e ©3^7 it is true: 

[337] = {1,1,1,1,1,1} 



w(3: 



37) 



^(337) = (0,0,4,2) 

[(3-0)!]° [(3-1)!]° _ 6" -2° 
1!1!1!1!1!1! ~ 1 
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For graphs 538 G ©3,7 and 539 £ ©3,7 we get: 

[538] = [539] = {1,1,2,2} 



'^(338) = ^(339) 

Then for the set ©3.7 we get: 



*(.938) = *(.939) - (0,1,2,3) 

[(3-0)!]° [(3-1)!]' 6" -21 1 



1! 1! 2! 2! 



22 



0(3,7)= J2 ^(5) = 1 + ^ + ^ = 2 



(21) 



3.2.8 fc= 



Graph 1740, which is displayed in Figure [T^ is the only bipartite graph belonging 
to the set ©3 § in the case n = 3 and fc = 8. 




Figure 12: n = 3, A: = 8 
For graph 540 G ©3,8 it is true: 

[ff4o] = {l,l,2,2} 

*(<74o) = (0,0,2,4) 
[(3-0)!]° [(3-1)!]" _ 60-20 _ 1 



w(54o) = 



1! 1! 2! 2! 



22 4 



Therefore: 



0(3,8)= ^(3) = I 



(22) 



3.2.9 fc = 9 



When n = 3 and fc = 9 there is only one graph and this is the complete bipartite 
graph (741 which is shown in Figure 1131 
For graph (741 is true: 

[541] = {3, 3} 
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Figure 13: n — 3, k — 9 



*(54i) = (0,0,0,6) 
[(3-0)!]° [(3-1)!]" 60-2" 



Therefore 



3! 3! 



0(3,9)= ""(a) 

sees, 9 



62 



1 

36 



1 

36 



(23) 



Having in mind the formula ([5]) and formulas (IT5|) -f- for the number Dg 
of all ordered pairs disjoint S-permutation matrices in n = 3 we finally get: 



Dg = (3!)i2 + (3!)« 



(24) 

1 1 
4 " 36 
= 1 260 085 248. 

The number dg of all non-ordered pairs disjoint matrices from Eg is equal to 



2 176 782 336+1 679 616 -1296 + 1008 - 352+ 125 - 37 + 8 -2 



dg = ^Dg = 630 042 624 



(25) 



3.3 On a combinatorial problem of graph theory related 
to the number of Sudoku matrices 

Problem 1 Let n > 2 is a natural number and let G be a simple graph having 
(n!)2" vertices. Let each vertex of G be identified with an element of the set E„2 
of all X S-permutation matrices. Two vertices are connected by an edge if 
and only if the corresponding matrices are disjoint. The problem is to find the 
number of all complete subgraphs of G having vertices: 

Note that the number of edges in graph G is equal to d„2 and can be calcu- 
lated using formula ([2]) and formula (jSj (respectively formulas (H)), ([71), ([8]) and 
©). 
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Denote by z„ the solution of the Problem [T] and let tT„ is the number of all 
r? X r? Sudoku matrices. Then according to Proposition [1] and the method of 
construction of the graph C?, it follows that the next equality is valid: 

^.-^ w 

We do not know a general formula for finding the number of all r? x r? 
Sudoku matrices for each natural number n > 2 and we consider that this is an 
open combinatorial problem. Only some special cases are known. For example 
in n = 2 it is known that (J2 = 288 ^7. . Then according to formula ([26| we get: 

^2 288 

In [4] it has been shown that in n = 3 there are exactly 

(73 = 6 670 903 752 021 072 936 960 = 

= 9! X 72^ X 2^ X 27 704 267 971 = 
2^" X 3^ X 5^ X 7^ X 27 704 267 971^ - 6.671 x 10^^ 
number of Sudoku matrices. Then according to formula we get: 



(73 6 670 903 752 021 072 936 960 „ ^oo 

Z3 = -f = ^TTHTTTTTT^ = 18 383 222 420 692 992 

9! 362 880 

References 

[1] R.A. Bailey, P.J. Cameron, and R. Connelly. Sudoku, gerechte designs, 
resolutions, affine space, spreads, reguli, and hamming codes. Amer. Math. 
Monthly, (115):383-404, 2008. 

[2] G. Dahl. Permutation matrices related to Sudoku. Linear Algebra and its 
Applications, (430):2457~2463, 2009. 

[3] R. Diestel. Graph Theory. Springer- Verlag Heidelberg, New York, 1997, 
2000, 2006. 

[4] B. Felgenhauer and F. Jarvis. Enumerating possible Sudoku grids. 2005. 
http: / / www.afjarvis.staff.shef.ac.uk/sudoku/sudoku.pdf. 

[5] R. Fontana. Fraction of permutations - an application to sudoku. Journal 
of Statistical Planning and Inference, (141):3697-3704, 2011. 

[6] F. Harary. Graph Theory. Addison- Wesley, Massachusetts, 1998. 

[7] H. Kostadinova and K. Yordzhev. An entertaining example for the usage of 
bitwise operations in programming. In FMNS-2011, volume 1, pages 159- 
168, Blagoevgrad, Bulgaria, 2011. SWU. 

[8] K. Yordzhev. On the number of disjoint pairs of s-permutation matrices. 
2012. arXiv:1211.1628. 



17 



